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ABSTRACT 


Multiple  coherence  gives  a  quantitative  measure  versus 
frequency  of  how  well  a  linear  combination  of  n  input  channels 
can  match  the  (n  +  1) st  channel  in  a  seismic  array.  If  the 
inputs  can  match  the  output  exactly,  then  the  multiple  coherence 
is  unity  and  only  n  channels  are  necessary  to  describe  the  noise 
field.  This  report  shows  multiple  coherence  versus  frequency 
with  2  to  9  input  channels  for  short  period  noise  fields  at  LASA. 

Intersubarray  noise  from  the  center  seismometers  in  the  500- 
foot  holes  at  LASA  shows  low  multiple  coherence  for  all  frequencies 
from  0.1  to  2.5  cps  even  with  the  closest  subarrays  represented. 

The  intersubarray  multiple  coherence  indicates  that  the  expected 
noise  reduction  from  a  prediction  error  filter  is  about  1  db 
over  the  fitting  interval  and  about  0  db  outside  the  fitting  in¬ 
terval. 

Intrasubarray  noise  has  a  high  multiple  coherence  for  fre¬ 
quencies  of  1  cps  and  lower  when  some  input  channels  are  within 
h  km  of  the  output  channel.  When  the  inputs  are  at  least  2  km 
away,  the  multiple  coherence  is  less  than  0.8  for  all  but  micro- 
seismic  frequencies. 


1  INTRODUCTION 


Most  basic  data  processing  techniques  for  signal  enhancement 
or  identification  depend  upon  the  structure  of  the  noise  within 
the  seismic  array .  If  some  of  the  coherent  noise  is  due  to  site 
characteristics  such  as  consistently  coherent  noises  from  par¬ 
ticular  directions,  then  techniques  using  multiple  coherence  .will 
help  to  isolate  these  consistent  linear  relations.  Many  optimum 
filters  for  estimating  the  signal  take  account  of  these  linear 
relations  implicitly  by  weighting  with  the  inverse  of  the  spectral 
noise  matrix.  However,  one  cannot  tell  whether  the  coherent  noise 
involved  is  due  to  noise  generating  events  which  cannot  be  pre¬ 
dicted  or  controlled.  Thus,  the  filters  must  be  recalculated  over 
a  period  of  noise  recording  immediately  prior  to  the  arrival  of 
each  single  signal.  Part  of  the  coherent  noise  generated  within 
the  array  may  be  due  to  various  causal  factors  for  c  particular 
array.  If  so,  we  can  learn  something  about  these  factors  by  ex¬ 
amining  the  linear  relations  between  the  various  array  elements. 

A  potential  benefit  here  is  that  a  consistent  linear  model  re¬ 
lating  the  different  sub-elements  would  eliminate  the  need  for 
computing  a  different  set  of  filter  coefficients  for  each  event. 

The  multiple  coherence  function  can  indicate  how  many  seis¬ 
mometer  outputs  in  an  array  are  necessary  to  properly  determine 
the  seismic  noise  field.  If  there  are  n  independent  seismic 
noise  components,  then  the  multiple  coherence  function  would  be 
unity  when  (n  +  l)st  seismometers  are  placed  in  an  array  to 
measure  seismic  noise  records.  If  part  of  the  background  is 
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composed  of  incoherent  noise,  then  the  multiple  coherence  function 
would  indicate  the  percentage  of  coherent  noise  present  and  the 
number  of  seismometers  necessary  to  define  this  coherent  noise. 

The  filter  relations  determined  by  the  multiple  coherence  com¬ 
putations  can  then  be  used  in  array  summation  to  bring  the  noise 

into  destructive  interference.* 

This  analysis  does  not  guarantee  that  such  optimum  processing 
is  possible.  For  example,  if  the  noise  and  signal  propagation 
characteristics  across  the  array  are  identical,  no  velocity  fil¬ 
tering  scheme  can  be  expected  to  separate  the  two  even  though 

the  multiple  coherence  might  be  unity. 

The  multiple  coherence  function  is  the  frequency  domain  equi¬ 
valent  of  the  prediction  error  filter  in  time.  If  n  input  seismic 
traces  predict  the  (n  +  1) st  trace  in  an  array  completely,  then 
the  multiple  coherence  will  be  unity  and  a  prediction  error  filter 
could  be  used  to  exactly  predict  this  (n  +  1) st  output.  In  fact, 
linear  filter  relations  derived  by  the  multiple  coherence  program 
produce  an  estimate  of  the  (n  +  1) st  trace  which,  when  subtracted 
from  the  actual  (n  +  1)  st  trace,  given  a  prediction  errot:  trace. 
Thus  the  combination  of  the  filter  derived  in  the  multiple  coher¬ 
ence  program  and  the  subtraction  operation  produces  a  prediction 
error  filter  as  shown  in  the  following  diagram. 


*For  the  mathematical  description  of  the  multiple  coherence 
computation,  see  Appendix  I. 
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The  first  objective  of  this  study  is  to  use  the  multiple 
coherence  function  to  estimate  the  degree  of  predictability  of 
the  short  period  noise  field  at  LASA.  The  result  in  turn  should 
tell  us  how  much  the  noise  power  should  be  reduced  by  optimum 
filtering  (e.g.,  maximum  likelihood  or  Wiener  filters)  if  the 
filters  were  theoretically  ideal. 

The  second  objective  is  to  determine  from  multiple  coherences 
the  number  of  independent  components  comprising  a  given  noise 
field  and  the  percentage  of  incoherent  noise  which  cannot  be  can¬ 
celled  by  any  kind  of  multichannel  filtering. 

The  third  objective  is  to  determine  the  stationarity  proper¬ 
ties  of  the  noise  field.  We  accomplish  this  by  applying  the 
»  multiple  coherence  program  to  three  different  time  samples  from 

the  same  array.  Then  the  multiple  coherency  filters  derived 


from  the  first  time  sample  are  applied  to  the  other  two  time 
samples.  If  the  filters  derived  in  the  first  time  sample  have 
done  a  good  job  of  predicting  the  noise  field  in  all  three 
samples,  then  the  data  are  said  to  be  stationary.  On  the  other 
hand,  if  the  filters  from  the  first  time  sample  do  a  progressively 
poorer  job  of  predicting  the  noise  in  the  other  samples  relative 
to  the  filters  associated  with  those  samples,  then  the  noise  is 
non-stationary  to  some  degree.  This  deterioration  in  predict¬ 
ability  of  the  multiple  coherence  filter  can  quantitively  measure 
the  non-stationarity  of  the  data**. 


**For  a  theoretical  discussion  of  the  stationarity  computation, 
see  Appendix  II. 
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2.  DESCRIPTION  OF  DATA 

We  computed  the  multiple  coherences  of  short-period  LASA 
data  for  both  within  the  subarray  (intrasubarray)  and  between 
subarrays  (intersubarray).  The  noise  within  a  subarray  was 
tested  for  subarray  AO.  The  output  seismometer  was  the  center 
element,  A010.  The  multiple  coherence  was  computed  as  a  function 
of  frequency  and  the  number  of  inputs.  The  first  inputs  added 
were  from  the  outermost  ring.  The  nunfcer  of  input  channels  in¬ 
creases  from  2  to  9.  The  multiple  coherence  is  plotted  versus 
frequency  with  the  number  of  input  channels  as  a  parameter. 

Since  adding  one  additional  channel  to  the  set  of  inputs  may  not 
increase  but  can  never  decrease  the  amount  of  input  information 
about  the  noise  field,  the  multiple  coherence  must  be  a  mono- 
tonically  increasing  function  with  increasing  number  of  inputs. 

In  computing  the  multiple  coherences  between  subarrays,  the 
center  elements  from  each  subarray  were  used.  Again  the  output 
trace  was  from  the  center  of  the  array,  subarray  AO.  The  order¬ 
ing  of  the  input  channels  was  done  in  two  ways.  In  the  first, v 
we  added  the  outer  most  channels  first  and  the  near  ones  last; 
in  the  second,  we  reversed  this  procedure. 

The  noises  examined  were  three  successive  signal-free  samples, 
eight  minutes  long,  recorded  on  25  March  1966.  In  addition, 
stationarity  tests  were  made  on  these  samples.  Also  tested  were 
one  signal-free  sample  and  one  sample  including  LONGSHCT  from 
29  October  1965. 
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3.  RESULTS 


Multiple  Coherences 

Figure  1  shows  intersubarray  multiple  coherences  versus 
frequency  for  two  time  samples  chosen  8  minutes  apart  for  25  March 
1966.  A  diagram  of  the  array  elements  chosen  is  shown  in  the  cen¬ 
ter  of  the  figure.  The  ordering  of  the  inputs  is  from  outermost 
to  closest  and  is  shown  in  the  figure.  The  multiple  coherences 
are  computed  every  .1  cps  over  a  range  from  .2  to  2.6  cps.  The 
number  of  points  in  each  sample  is  4799.  The  number  of  lags  com¬ 
puted  in  the  correlation  function  is  200.  The  data  sampling  rate 
is  10  points  per  second. 

The  intersubarray  multiple  coherences  are  seen  to  be  low 
over  the  entire  frequency  range  for  even  as  many  as  9  inputs. 

The  number  of  lags  used  in  the  correlation  is  sufficient  to  detect 
any  propagating  noise  component  with  a  velocity  greater  than  1.25  km 
per  second.  This  velocity  limitation  may  affect  the  correlation 
between  AO  and  the  D-ring.  However,  most  of  the  inputs  ore  closer 
than  15  km  so  that  the  lack  of  sufficient  lags  to  detect  signifi¬ 
cant  correlations  is  not  a  cause  for  the  low  multiple  coherences 
observed. 

Figure  2  shows  the  multiple  coherences  versus  frequency 
for  both  inter-  and  intra subarray  noise  from  a  sample  on  29  October 
1965.  The  subarray  and  array  diagrams  are  shown  in  the  center  of 
the  figure,  and  the  ordering  of  the  Inputs  is  again  from  the  outer¬ 
most  to  the  closest. 

The  intersubarray  noise  still  shows  low  coherences  although 
the  values  for  all  10  inputs  are  somewhat  above  those  plotted  in 
Figure  1.  The  intra subarray  coherences  are  quite  large  out  to  1  cps 
when  7  or  more  inputs  are  used.  However,  the  7th,  8th  and  9th 
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inputs  are  only  h  km  away  from  the  output  trace  so  high  coherences 
are  to  be  expected .  When  the  closest  inputs  are  2  km  away,  however, 
the  multiple  coherences  are  significantly  less  than  .8  for  almost 
all  frequencies.  Here  again,  the  number  of  lags  is  sufficient  to 
detect  any  propagating  noise  components  that  could  exist. 

Figure  3  shows  the  multiple  coherences  for  the  same  inter¬ 
subarray  data  given  in  Figure  2.  In  this  case,  the  ordering  of  the 
inputs  adds  the  closest  seismometers  first.  The  final  multiple 
coherence  trace  with  all  10  inputs  must  still  be,  and  is,  the  same 
as  the  final  trace  on  Figure  2.  We  find,  however,  that  even  using 
the  closest  seismometers  first,  we  do  not  approach  this  final  mul¬ 
tiple  coherence  any  faster.  This  result  shows  that  the  intersub¬ 
array  noise  at  LASA  is  essentially  incoherent. 

As  a  demonstration  that  the  multiple  coherence  program 
will  detect  significant  correlations  when  they  do  exist,  we  have 
shown  the  intrasubarray  multiple  coherence  when  the  noise  field 
contains  a  seismic  signal.  The  upper  diagram  in  Figure  4  shows 
the  multiple  coherence  within  a  subarray  when  the  LONGSHOT  event 
was  recorded.  Here  the  multiple  coherence  is  above  .8  for  most 
frequencies  with  5—6  inputs,  and  is  essentially  unity  with  9 
inputs.  The  lower  diagram  in  Figure  4  shows  similar  multiple 
coherences  when  the  data  contains  a  smaller  seismic  signal.  Here 
the  multiple  coherence  is  reduced  because  the  signal  is  weaker 
compared  to  the  incoherent  background  than  in  the  case  of  LONGSHOT. 
Even  so,  the  multiple  coherence  for  this  weaker  signal  plus  noise 
is  significantly  above  the  multiple  coherence  for  noise  alone  with¬ 


in  a  subarray. 


Figure  5  shows  intersubarray  multiple  coherences  for 
the  LONGSHOT  and  Andreanof  events  similar  to  the  intrasubarray 
coherences  of  Figure  4.  The  multiple  coherences  on  Figure  5 
show  striking  differences  from  those  of  Figure  4.  The  intra¬ 
subarray  coherences  are  highest  at  the  low  frequencies  and  some¬ 
what  lower  over  the  middle  frequencies.  On  Figure  5  the  lowest 
frequencies  show  low  coherences  and  quite  high  coherences  around 
the  middle  frequencies.  This  result  can  be  esqolained  by  con¬ 
sidering  the  signal-to-noise  ratios.  At  the  lowest  frequencies 
we  have  noise  only.  Over  intrasubarray  distances  the  microse..smic 
noise  is  coherent  whereas  over  intersubarray  distances  it  is  not. 
Over  the  middle  frequencies  where  the  signal  spectra  are  strongest, 
the  noise  is  weakest  for  the  intersubarray  case  since  the  inputs 
are  all  from  500  foot  holes.  Consequently,  signal-to-noise  ratios 
from  the  intersubarray  case  especially  for  the  Andreanof  event, 
is  much  higher  and  the  multiple  coherences  consequently,  much 
higher  than  for  the  intra subarray  events. 

Power  Spectra 

Spectra  for  the  multiple  coherence  examples  on  Figures  1-4 
are  presented  on  Figures  6-10  and  Figures  12  and  13.  Figure  6 
shows  the  output  spectra  at  AO  and  the  range  of  input  spectra  for  9 
intersubarray  inputs  from  the  first  noise  sample  on  March  25,  1966. 

Figure  7  snows  the  spectra  from  the  same  500-foot  holes 
shown  in  Figure  6  but  computed  with  100  instead  of  200  lags.  The 
variation  in  lags  shows  essentially  no  difference  in  the  power 
spectra. 

In  Figure  8  the  power  spectra  shown  for  the  third  time 
sample  are  similar  to  those  shown  for  the  first  in  Figure  6. 
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Figure  9  shows  the  output  and  input  power  spectra  for  the 
intrasubarray  noise  sample  immediately  preceding  the  LONGSHOT 
event.  Figure  10  shows  the  similar  spectra  for  the  intex subarray 
noise  prior  to  LONGSHOT.  The  two  output  traces  in  Figure  9  and  10 
are  identical. 

A  marked  difference  between  the  input  spectra  for  the 
higher  frequencies  shows  up  between  Figures  9  and  10.  The  intra¬ 
subarray  spectra  (Figure  9)  are  from  200-foot  holes  and  show 
appreciably  higher  energy  near  2  cps  than  do  the  intersubarray 
spectra  from  500— foot  holes  shown  in  Figure  10 .  The  spectral  dif¬ 
ference  between  the  20C  and  500  foot  seismograms  implies  that  the 
intrasubarray  multiple  coherence  might  be  improved/  particularly 
at  the  higher  frequencies,  if  we  would  use  an  output  trace  from 
a  200-foot  hole.  Figure  11  shows  the  intrasubarray  multiple  co¬ 
herence  when  the  output  is  from  a  200-foot  hole  h  km  from  the 
center  seismometer.  Except  for  this  difference  in  output  traces 
this  case  corresponds  to  that  shown  in  the  upper  half  of  Figure  2. 
In  comparing  the  two  we  find  virtually  no  difference  whatever 
between  the  multiple  coherences  for  any  frequency  in  these  two 
cases.  Furthermore,  the  multiple  coherences  versus  frequency 
for  these  two  cases  agree  regardless  of  the  number  of  input 
channels  used. 

Figure  12  shows  the  output  and  input  spectra  for  the 
interval  including  the  LONGSHOT  event. 

Figure  13  shows  the  input  and  output  spectra  for  the 
interval  including  the  Andreanof  earthquake. 
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Stationarity  Tests 

The  multiple  coherence  program  derives  a  set  of  n  filters 
for  the  n  input  seismograms  which  together  provide  the  best  linear 
estimate  for  the  (n  +  1) st  seismic  trace.  The  difference  between 
the  observed  (n  +  l)st  trace  and  the  best  estimate- is  the  error 
trace,  if  the  multiple  coherence  is  unity,  the  prediction  is  per¬ 
fect  and  the  error  trace  will  be  zero.  If  we  form  the  ratio  of 
the  error  spectra  over  the  observed  spectra,  we  can  get  a  measure 
of  the  reduction  in  noise  power  possible  from  the  theoretical 
optimum  filters.  Thus  the  db  improvement  as  a  function  of  fre¬ 
quency  can  be  expressed  as 

db  =  10  log  (error/observed) . 

The  prediction  error  filters  will  do  the  best  job  in 
eliminating  the  noise  background  when  they  are  applied  to  the 
noise  sample  from  which  they  are  derived.  However,  if  the  noise 
is  stationary,  the  same  filter  could  be  expected  to  do  nearly  as 
well  when  applied  to  later  time  samples  from  the  same  array. 

Figure  14  shows  the  expected  noise  reduction  in  db  when  the  pre¬ 
diction  error  filters  that  were  derived  from  the  first  time  sample 
are  applied  to  the  first  time  sample.  In  addition  these  same 
filters  are  applied  to  the  second  and  third  time  samples.  The 
three  8-minute  time  samples  are  adjacent  to  each  other.  Figure  14 
shows  that  the  reduction  of  the  noise  background  from  a  prediction 
error  filter  can  be  expected  to  be  between  .5  and  2  db.  The  same 
filter  applied  to  the  second  and  third  time  sample  gives  0  db  and 
thus  shows  no  reduction  of  the  noise  background  whatever.  This 
small  amount  of  noise  reduction  of  intersubarray  data  is  to  be 
expects Ct  when  the  multiple  coherence  is  less  than  30%. 


-10- 


Figure  15  shows  similar  results  when  filters  derived 
from  the  second  time  sample  are  applied  to  all  three.  This 
expected  noise  reduction  and  lack  of  stationarity  is  further 
evidence  that  the  noise  is  essentially  incoherent. 
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4.  CONCLUSIONS 


1.  Inter subarray  noise  at  LASA  shows  a  low  multiple  coherence 
over  all  frequencies,  even  with  .9  input  channels. 

2.  Even  che  closest  subarrays  at  LASA  show  little  coherency. 

3.  Intrasubarray  noise  at  LASA  shows  a  high  multiple  coher¬ 
ence  for  frequencies  of  1  cps  and  lower  when  7  or  more  inputs  are 
used.  However,  the  7th,  8th,  and  9th  inputs  are  only  h  km  away 
from  the  output  trace.  When  the  closest  inputs  are  2  km  away  from 
the  output,  the  multiple  coherences  are  significantly  less  than  .8 
for  all  but  microseismic  frequencies, 

4.  The  digital  programs  will  show  high  coherence  and  multiple 
coherence  when  present  e.g.,  LONGSHOT  event  has  high  multiple  co¬ 
herence  over  subarrays.  An  Andreanof  earthquake,  a  weaker  event 
than  LONGSHOT,  shows  multiple  coherences  between  those  of  LONGSHOT 
and  noise-only  cases. 

5.  Low  multiple  coherence  is  not  due  to  insufficient  lags. 

The  programs  have  sufficient  resolution  to  handle  noise  propagation 
with  velocity,  1  km/sec  and  less. 

6.  Power  spectra  from  the  seismometers  at  500  feet  agree  with 
those  at  200  feet  only  out  to  1.2  cps.  Beyond  this  frequency  the 
power  spectra  density  of  the  deeper  seismometers  continue  to  fall, 
while  the  power  spectra  densities  from  the  shallow  seismometers 
rise. 

7.  The  multiple  coherence  using  shallow  seismometers  as 
inputs  trying  to  match  the  deeper  seismometer  as  output  falls 
rapidly  beyond  1.2  cps. 

8.  The  intra subarray  multiple  coherence  using  shallow  seis¬ 
mometers  as  both  inputs  and  output  is  approximately  the  same  as 
for  intrasubarray  multiple  coherences  which  uses  the  trace  from 
the  500-foot  hole  as  output.  This  result  is  obtained  in  spite 
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of  a  striking  difference  in  the  power  spectra  of  the  shallow  and 
deep  seismograms. 

9.  The  inter subarray  multiple  coherence  indicates  that  the 
expected  noise  reduction  from  a  prediction  error  filter  is  about 
1  db  over  the  fitting  interval  and  about  0  db  outside  the  fitting 
interval. 
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Figure  2.  Multiple  coherences  versus  frequency  for  both 
inter  and  intrasubarray  noise  recorded  at  LASA 
on  29  October  1965. 
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Figure  4.  Intrasubarray  multiple  coherences  versus  frequency 
for  the  LONGSHOT  event  and  an  Andreanof  earthquake 
event . 
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Figure  5.  Intersubarray  multiple  coherences  versus  frequency 
for  the  LONGSHOT  events. 
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Figure  6.  Output  and  range  of  input  spectra  for  intersubarray 
samples  used  on  Figure  1. 
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Figure  9.  Output  and  range  of  input  spectra  for  the  intra¬ 
subarray  noise  immediately  preceding  the  LONGSHOT 
event . 
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Figure  10. 


Output  and  range  of  input  spectra  for  the  inter¬ 
subarray  noise  immediately  preceding  the  L0NGSH0T 

event . 
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Figure  11. 


Intrasubarray  multiple  coherence  versus  frequency 
when  the  output  is  from  a  seismometer  in  a  200 
foot  hole  instead  of  the  500  foot  hole  at  the  sub¬ 
array  center. 
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Figure  13.  Output  and  range  of  input  spectra  for  the  interval 
including  the  Andreanof  earthquake  event. 
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appendix  I 

*Multiple  Coherence  Functions 


Consider  a  collection  of  q  clearly  defined  inputs  x^t); 

i  -  1,2,...  ,q,  and  one  output  y(t),  as  pictured  in  Figure  5.12. 
Let  ©i  (f)  =  Gii(f)  be  the 


Figure  5.12  Multiple-input  linear  system. 

power  spectral  density  function  for  x^t),  and  Gi;J<f)  be  the 
cross-spectral  density  function  between  x^(t)  and  x^(t).  De¬ 
fine  the  N  x  N  spectral  matrix  by 


Gn( f)  Gu(f)  •••  Gu(f)~ 
c'n(f)  Gn(f)  Gh(f) 


L^iC 0  GJJ)  GJif) J 


(l) 


^  rhis  explanaticn  of  multiple  coherence  functions  was  taken  from 
Measurement  and  Analysis  of  Random  Data",  Bendat ,  J.  s.,  and 
Piersol ,  A.  G.,  John  Wiley  end  Sons,  1966.  For  more  detailed  theo- 

partui  “a 


i-i 


The  ordinary  coherence  function  between  x^(f)  and  x^(t)  is 
defined  by 


G&f)  Gt(f) 


(2) 


The  multiple  coherence  function  between  x^ft)  and  all  /ther 
inputs  x1(t),  x2  (t) .excluding  x^t),  is  defined  by 

v,  AD- 1  - \OAJ)  c‘(/)r‘  ,3) 

where  G1(g)  denotes  the  ith  diagonal  element  of  the  inverse 
matrix  G^ff)  ^associated  with  Eq.  (1).  The  ordinary  and 
multiple  coherence  functions  are  both  real-valued  quantities 
which  are  bounded  by  zero  and  unity.  That  is, 


o  £  r<W  £  i  (4) 

o  <;  yLU  )  <.  i 

The  multiple  coherence  function  is  a  measure  of  the  linear 
relationship  between  the  time  history  at  one  point,  and  the  time 
histories  at  the  collection  of  other  points.  That  is,  the  mul¬ 
tiple  coherence  function  indicates  whether  or  not  the  data  at 
all  of  the  other  points  linearly  produce  the  results  at  a  given 
point . 
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APPENDIX  2 


Theoretical  Development  of  The  Stationarity  Relations 

A.  Noise  Reduction  Within  The  Fitting  Interval 

A  number  of  useful  statistical  measures  such  as  ordinary 
and  multiple  coherence  can  be  used  as  tools  to  indicate  the  amount 
of  noise  reduction  feasible  in  a  multiply  coherent  array.  The 
basic  linear  model  vdiich  determines  the  db  reduction  possible  in 
the  noise  field  by  multiple  coherence  filtering  relates  a  refer¬ 
ence  element  (trace)  y(t)  of  an  array  to  the  other  elements,  say 
x^t),  x2(fc),  xp(t)  in  the  array  through  the  linear  model 


y(t)  - 


P 


L 


00 

J 


h^ (nr.)  xk  (t -nr.)  da 


(1) 


k-1 

Generally  we  determine  h^t)  as  \:he  time  invariant  linear  filter 
that  makes  the  mean  square  error  between  y(t)  and  its  predicted 
value  a  minimum,  i.  e. 


P  - 

E  |  y(t)  -  £  J  hk  (a)  xk(t-a)  da|  -  min  (2) 

k-1  -• 

which,  by  the  usual  orthogonality  principle  (see  Papoulis  1), 
yields  the  condition 

P 

Ey(t)  x^(t+T)  -  £  J  h^a)  E  Xj^  (t-a)  x^  (t+T)  da 
k-1 

l  -  1,2,  ...,  P,--<T<»  (3) 


2-1  - 


♦ 


or 

v(t)  •£  J  va>  V*. <w0 


k-l 


Yt 


da 


(4) 


Which  by  taking  Fourier  transforms  implies  that 
P 

(U)) 


(5) 


V*1"’  "A  \*t 

k-l 

Ncv,  the  mean  square  error  can  be  written 

P  p  «  N 

E|y(t)  -  I  J  lya)  xk(t-a)da|2  -  E  (y(t)  -  £  J  h^ajx^t-ejdaj  y(t) 


k-l 


k-l  — 


Ryy(0)  “  ^  J  Va)  \y{tl)  dft 


k-l 


-  J  [  V’  -  I  "*  (“”  V(“’  ^ 


duj 


2tt 


k-l 


KS  -  8  S_1  S  ")  (ID)  — : 

yy  ^y*  xx  ~xy  / 


das 

2tt 


.1(1-  «2  <»>  )  Syy(tt,) 


du) 


2tt 


(6) 


Where  a2 (id)  is  the  multiple  coherence  and  (1  -  a  (id)  )  measures 

2 

the  reduction  in  power  possible  at  the  frequency  to.  With  a  (u>)  -  1 

2 

the  mean  square  error  is  zero  and  with  a  (u>)  ■  0  the  mean  square 
error  is  just 


J  V*’  "5^  "  Ryy  <0)  ’  E  ly 


(7) 


vtfiich  is  the  original  power  in  the  process  y(t) .  Now  the  quantity 
(l-a^(u))  )  Syy(tt>)  represents  the  residual  power  at  each  frequency 
after  the  best  linear  estimate  of  che  form  (1)  has  been  subtracted 
out.  Hence  the  db  reduction  in  power  at  each  frequency  is  jt.st 
the  ratio  of  the  output  power  of  the  residual  (see  equation  (2)  ) 
to  the  input  power  in  y(t)  or 


See  {m)  2 

IH (u>)  -  10  log  -g - ^ —  *  10  log  (1  -  a  (u>)  ) 


(8) 


2 

where  n  Im)  is  the  multiple  coherence  and  S  (tv)  is  the  power 

ee 

spectrum  of  the  error  process 
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y(t)  -  I  J  h^(a)  xk(t  -  or.)  da 


k-1  -« 


Noise  Reduction  Outside  The  Fitting  Interval 

We  would  also  like  to  determine  the  noise  reduction  in  db 


tAiich  would  result  from  using  a  set  of  filters  g^(t),  k"l,  . . .  ,p 
which  have  been  derived  either  from  another  fitting  interval  or 


from  theoretical  considerations.  To  accomplish  this  let  h^tt) 
be  the  optimal  filters  for  the  time  under  investigation  and  let 


g^(t)  be  any  other  set  of  filters  whose  mean  square  error  is  to 
be  compared  with  h^(t).  The  mean  square  error  of  the  g  filters 
can  be  written  using  the  orthogonality  principle  as 
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Hence,  if  we  call  the  new  error  t‘ (t)  we  have 


P 

r>  r 


V  (t)  =  Y(0  -  £  J  gk(a)  x  (t-a)  da 


k=l 


with  power  spectrum  Sfl/ ,(«,),  the  new  value  for  the  improvement 


in  the  sk(t)  filters  would  be 
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yy 
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Equation  (12)  shows  that  the  improvement  in  the  3k(t)  filters  is 
expressed  in  terms  of  the  improvement  in  the  hk(t)  filters  and  a 
correction  term  which  is  zero  vdien  H  =  G. 


(12) 
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The  improvement  values  yiu)  and  IG(u»  in  equations  (8)  and  (12) 
are  those  shown  in  the  main  body  of  the  report. 
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